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ABSTRACT 

The dynamical equations describing the evolution of a physical system generally have a free- 
dom in the choice of units, where different choices correspond to different physical systems 
that are described by the same equations. Since there are three basic physical units, of mass, 
length and time, there are up to three free parameters in such a rescaling of the units, Nf 4.3. 
In Newtonian hydrodynamics, e.g., there are indeed usually three free parameters, Nf = 3. 
If, however, the dynamical equations contain a universal dimensional constant, such as the 
speed of light in vacuum c or the gravitational constant G, then the requirement that its value 
remains the same imposes a constraint on the rescaling, which reduces its number of free 
parameters by one, to Nf = 2. This is the case, for example, in magneto-hydrodynamics or 
special-relativistic hydrodynamics, where c appears in the dynamical equations and forces the 
length and time units to scale by the same factor, or in Newtonian gravity where the gravita- 
tional constant G appears in the equations. More generally, when there are A^ u( jc independent 
(in terms of their units) universal dimensional constants, then the number of free parameters 
is Nf - max(0, 3 - N u d c ). When both gravity and relativity are included, there is only one free 
parameter (Nf — 1, as both G and c appear in the equations so that N u a c = 2), and the units 
of mass, length and time must all scale by the same factor. The explicit rescalings for dif- 
ferent types of systems are discussed and summarized here. Such rescalings of the units also 
hold for discrete particles, e.g. in A^-body or particle in cell simulations. They are very useful 
when numerically investigating a large parameter space or when attempting to fit particular 
experimental results, by significantly reducing the required number of simulations. 

Key words: methods: numerical — methods: miscellaneous — methods: N-body simulations 
— hydrodynamics — magnetohydrodynamics (MHD) — gravitation 



1 INTRODUCTION 

Numerical simulations are gradually but steadily playing an in- 
creasingly more important role in the study of many physical sys- 
tems. In particular, they have a growing impact on many different 
fields, such as fluid dynamics, plasma physics, astrophysics, par- 
ticle physics, Earth and planetary sciences, and meteorology. In 
particular, numerical computation is a vital tool for studying com- 
plex problems that are hard to solve analytically, such as nonlinear, 
many-body or multiple scale processes. 

Here I outline how for many types of numerical simulations, 
the results of a single simulation correspond to a whole family of 
physical systems. This arises from the more general freedom in the 
choice of units in the dynamical equations that describe a particular 
type of physical system. This freedom (or lack thereof) also holds 
for solutions of these equations, as long as the initial or bound- 
ary conditions do not impose additional universal dimensional con- 
stants (UDCs), thus increasing N lK \ c . It is equally valid for analytic 



or numerical solutions of these equations. Here the focus is on nu- 
merical solutions, and in particular on numerical simulations. 

In this work it is assumed that the number of basic physical 
units is iVbpu = 3, corresponding to units of mass (m), length (/) and 
time (f). This means that electric charge is expressed in terms of 
these units ([q] = m l/2 P /2 t~ l ), which corresponds to c.g.s Gaussian 
units in which Maxwell's equations contain one UDC - the speed of 
light in vacuum, c. The number of free parameters Nf that span the 
family of physical systems that correspond to a single simulation is 
generally given by 



Nf = max (0, N bpu - N udc ) 



(1) 



For magneto-hydrodynamics (MHD), either Newtonian or special 
relativistic, where the only UDC is c (i.e. N udc = 1), this implies 
that Nf = 3-1=2. 

An alternative choice of basic physical units is SI units, in 
which electric charge is explicitly treated as a fourth basic physi- 



cal unit (measured in Coulombs), so that N\ 



bpu 



4. In these units, 
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however, Maxwell's equations contain two UDCs, e (the permit- 
tivity of free space or the electric constant) and yUo (the permeabil- 
ity of free space or the magnetic constant), where (e p )~ l/2 = c - 
Therefore, with this choice of units _/V u j c = 2 for MHD, leading 
to Nf = 4 - 2 = 2, i.e. the same number of free parameters as 
before. This nicely demonstrates that the freedom in rescaling the 
basic physical units for a given physical system is independent of 
the specific choice of basic units. Therefore, the particular choice 
that we use in this work should not affect Nf and the implied rela- 
tions between the allowed scalings of the units of mass, length and 
time. 

The concept of UDCs is not new. For example, lEllisI Jl968h 
has referred to UDCs as universal scale-dependent constants, and 
demonstrated that if units of mass and force are expressed in terms 
of length and time ([;;?] — » l 3 f~ 2 , [F] — » / 4 r~ 4 ) then this renders 
the gravitational constant G dimensionless in such units, so that an 
appropriate choice of their magnitude can make it equal to unity 
and thus disappear from Newton's law of gravity. In such units, 
Mick = 2 (as the basic physical units are only of length and time) 
while iV„d C = for Newtonian gravity, leading to Nf = 2 - = 2. 
For the regular units A'bpu = 3 and N mk = 1 for Newtonian grav- 
ity, again leading to Nf = 3 — 1 = 2. Unlike the c.g.s Gaussian 
units that were constructed in order to eliminate the constant [or 
SI UDC l/(4ne )] in Coulomb's law, the above choice of units that 
eliminate G from Newton's law of gravity do not work very well in 
general. The reason for this is the equivalence (which does not have 
an electromagnetic analog) between gravitational mass, for which 
these units were constructed, and inertial mass, which appears also 
in systems in which gravity is unimportant and can be neglected. In 
such systems that choice of units will artificially eliminate a degree 
of freedom in the rescaling of the units for no good reason. 

In §[2] the scalings are explicitly derived for different types of 
simulations, and summarized in Table[T] Some caveats, namely the 
possible increase in N uic for certain equations of state or radiative 
processes, are discussed in § [3] In § |4] the application of special 
relativistic hydrodynamic simulations to model gamma-ray burst 
(GRB) afterglow is discussed in some more detail, as a useful case 
study. The conclusions are discussed in §|5] 



2 SCALING RELATIONS FOR DIFFERENT TYPES OF 
NUMERICAL SIMULATIONS 

Physical systems are usually described either in the contin- 
uum limit, such as in hydrodynamics or magneto-hydrodynamics 
(MHD), or by following the motions of discrete particles. Numer- 
ically, the former description corresponds to hydrodyna mic (e.g., 
iLi et al.ll2O10l; lLawson & Barakosll2oTll: IWallace 201 H) or MHD 



(e.g. jFendt & Memold2008l : llshihara et alj2009l : lAmit et alj20ld) 

simulations, while the latter includes examples such as cosmo- 
logical A'-body simulations (with a large number N » 1 point 
particles that interact only through their mutual gravitational at- 



traction; 



eg 



Kravtsov et aill 19971: iNavarro. Frenk & White! 1 19971 ; 



ISpringeletal.ll2005l) or particle in cell (PIC) plasma simulations 
(where a large number of point particles with both positi ve and neg- 
ative electric charges interact electromagnetically; e.g. Birn et all 
1200 ll : IPukhov & Mever-ter-VehrJl2002l ; ISpitkovskvll2008l) . Table [Tl 
summarizes the allowed scalings for different types of simulations. 
Hybrid simulations that include both a continuous medium and dis- 
crete particles are also possible (e.g.. | Katz. Weinberg & Hernquisll 
1 19961 : iGnedin et al.ll2004l : ISpringelll2005h . and the restrictions on 
them can be derived in a straightforward manner by combining 



the restrictions on their constituents. This is manifested in the total 
number of independent (in terms of their units) UDCs, N^, which 
determines the number of free parameters Nf that span the family 
of physical systems that correspond to a single simulation through 



N f = max(0, 3-jV udc ) , 



(2) 



where we assume for the rest of this work that there are three basic 
physical units (A'bpu = 3) of mass, length and time. 

Hydrodynamic Simulations: here we outline how the results of 
different types of hydrodynamic simulations correspond to a fam- 
ily of physical systems. The numerical code solves, e.g., for the 
proper rest-mass density p, pressure p, and velocity v = [3c, as a 
function of time and space, (t, r), and assumes some equation of 
state that relates the specific enthalpy to the pressure and density. 
The evolution of these quantities is usually solved over a finite vol- 
ume V and time range /, < r < tf. The initial conditions at r, must 
be specified over the volume V, as well as boundary conditions at 
the edges of this volume at f,- < t < tf. In order to solve the hydro- 
dynamic equations numerically, they are first made dimensionless 
by moving to code units that are determined by choosing some spe- 
cific scales (which we shall denote by a subscript "0") for the three 
basic physical units of mass (ra ), length (/ ), and time (t ). The 
corresponding dimensionless variables in code units are denoted 
by a twiddle, where a general quantity Q with units of m A J B t c cor- 
responds to Q = Q/(m A l B t^). 

Once a particular initial physical configuration is mapped onto 
the dimensionless code variables, the hydrodynamic equations are 
solved for these variables. Then, the numerical solution is usually 
translated back to the original physical units, Q = Qmil B t9. This 
is not a unique procedure, however, because of the freedom in the 
choice of units that was described above. Therefore, the same nu- 
merical solution also holds equally well for a whole set or family 
of different physical systems, which correspond to different choices 
for the basic physical units, (m' , 1' , f ) = (fmo, &lo, J/'o)- Such 
different choices of units can conveniently be implemented when 
switching back from the dimensionless code variables to the corre- 
sponding variables with physical units. 

For a purely Newtonian simulation without gravity, there are 
indeed three free parameters (£, a and 77), i.e. Nf = 3 since there are 
no relevant UDCs (7Y u( jc = 0), and the family of physical systems 
represented by the simulation is given by 

Q' = Q(m' ) A (l' ) B (t' f = S A a B i 1 c Qm A l B t C Q = f*«V0 

<^> Q' (t' = r)t, r' = af) = ( A a B rj c Q(t, f) . (3) 

The scaling r' = af reads for different coordinate systems, 

(x , y', z ) = (ax, ay, az) cartezian , 
(z', r' j, ff) = (az, axqri, 0) cylindrical , 
(/, 6', (f)') = (ar, 9, <p) spherical . (4) 

Similarly, the initial conditions would be at (t' ,r') = (r]tj,ar) and 
the boundary conditions would be at the edge of V = a 3 V or (t\ < 
f < f f , r* cdge ) where r* cdgc = af cigs ,, t\ = T]tj and t' f = rjt f . 

When there are relativistic velocities, either of bulk motions or 
random motions of the particles (i.e. relativistic temperatures), then 
we no longer have/3 <K 1 and pip <k c 2 , so that relativistic effects in 
the dynamical equations (which depend on the bulk Lorentz factor 
T) or the equation of state (which depends on pip) can no longer be 
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neglected. This requires thaflr' = r (i.e. v' = v) and p' lp' = pip, 
where unprimed quantities are the usual ones for £ = a = t] = 1 , 
which implies that 77 = a, i.e. that the length and time units scale by 
the same factor. A more elegant way of deriving this is that special 
relativity introduces c (the speed of light in vacuum), a UDC with 
units of lit, and thus requires a/rj = 1 in order to keep its value 
the same in all our family of physical systems. This reduces the 
number of free scaling parameters to two (Nf = 2 since N u a c = 1). 
In particular, x'^ = (f,r) = a(t,r) = ax^p'jp = p'jp = (a~ 3 , 
and the family of physical systems corresponding to a particular 
simulation is given by 

K, 1' , t' ) = (Om, ah, otk) , Q'(x? M = ax,) = ? a B+c Q(x„) . (5) 

In particular, the total energy E, and mass M either in a particular 
computational cell or in the whole computational box (or volume 
V) scale as E'/E = M'/M= £. 

When gravity is included, this introduces the gravitational 
constant G - a UDC with units of m~ [ l 3 r 2 . Thus, in order for it 
to keep the same value in all our family of systems requires that 

(OVrVr 2 = Pm-'r 1 ^=> <r = cry 2 - (6) 

When there is only weak or Newtonian gravity (and general 
relativistic effects can be neglected), and Newtonian (bulk or ther- 
mal) motions, then G is the only UDC (A'udc = 1, Nf = 2) and 

(m' , 1' , t' ) = (a 3 r]~ 2 m , al , rjt ) , 
Q (t' = r\t, r> = ar) = a B+3A rf- 2A Q(t, r) . (7) 

When there effects of general relativity cannot be ignored, or 
for Newtonian gravity with relativistic velocities (either bulk or 
thermal), then there are two relevant UDCs, G and c (N uic = 2), 
which imply Eq. ® and a = t], respectively. Together this implies 
that f = a = T], i.e. that all three scaling coefficients are equal, 

(m' , 1' , t' ) = a (m , l , t ) , Q'(x[, = ax,) = a A+B+c Q(x„) . (8) 

Magneto-Hydrodynamic (MHD) Simulations: the MHD equa- 
tions are also based on Maxwell's equations, and thus include c as 
a UDC, so they require that a = r\. This holds even in the Newto- 
nian case, where there are two free parameters (Nf = 2) describing 
the relevant family of physical systems corresponding to a particu- 
lar simulation, according to Eq. 10. If gravity is included, even if 
weak or Newtonian gravity, then this introduces a second UDC, G, 
resulting in only one free parameter describing the relevant family 
of physical systems (Nf = 1 since A'udc = 2), according to Eq. 0. 

Simulations with discrete particles: there are various types of 
simulations that aim to describe the motions of discrete point-like 
particles, under the influence of the mutual forces that they exert 
on each other, rather than a continuous medium that is described 
by hydrodynamic or MHD equations. Here I briefly go over two 
important types of such simulations. 

The first type is particle in cell (PIC) simulations of the mo- 
tions of charged particles of either positive or negative electric 
charge under the mutual electromagnetic forces that they exert on 
each other. In this case Maxwell's equations introduce c as a UDC 
(implying a = rj). If we do not mind that the scaling would change 
the rest mass and/or electric charge of particles, then this would be 

1 For r» 1 there exists a different type of res caling that allows r* to vary 
relative to T I Mimica, Giannios & Aloy 2009), but this is not a rescaling 
of the basic physical units, and it is valid only for a forward-reverse shock 
system in which the forward shock is ultra-relativistic. 



Table 1. The freedom in the choice of the basic physical units for different 
types of numerical simulations. The columns, from left to right, list the type 
of simulation, the relevant independent (in terms of their units) universal 
dimensional constants (UDCs), their number A uac , the number of free pa- 
rameters Nf they allow for a rescaling of the units, and the imposed relation 
between the rescaling factors for mass (f ), length (a) and time (;;) units. 
' If one or more of the cosmological parameters (such as Ho or o"s) are 
treated as UDCs this reduces Nf - see the discussion in §fj] 



type of simulation 


UDCs 


A'udc 


N f 


constraints 
on rescaling 


Newtonian hydrodynamics 







3 




relativistic hydrodynamics; 
Newtonian/relativistic MHD 


c 


1 


2 


a = rj 


Newtonian gravity (e.g. in 










stellar/planetary dynamics, 


G 


1 


2 


c = «V 2 


cosmological A-body ' ) 










general relativistic 










hydrodynamics or MHD; 
Newtonian gravity + MHD 


G,c 


2 


1 


£ = a = r] 


or relativistic velocities 










particle in cell (PIC); 


c 


1 


2 


a = rj 


PIC + particular particles 


c, q, m 


3 





f = « = !?=! 



the only constraint (A'udc = 1), implying Nf = 2 and Eq. lf5J. If, 
however, it is important for us to accurately model a specific par- 
ticle species (such as electrons/positrons) of a given universal rest 
mass and electric charge, then this would add two more constraints 
(and altogether A/ U d c = 3), thus removing all the remaining freedom 
in the scaling parameters (Nf = 0) and implying f = a = rj = 1. 

The second type is A'-body simulations, that are often used in 
cosmology and stellar or planetary dynamics, where N point-like 
masses move under their mutual gravitational forces. Since gravity 
is the only force involved, G must obviously remain constant, im- 
plying Eq. {6]l. If there are only Newtonian gravity and velocities 
then there are two free parameters (Nf = 2, N ulic = 1) and Eq. 10 
holds. Otherwise, if relativistic effects cannot be neglected, then c 
also enters the relevant equations as a second UDC (A'udc = 2) re- 
sulting in only one free parameter (Nf =1), and implying Eq. ([SJ. 

In cosmological A'-body simulations with Newtonian gravity 
and velocities, the only bona fide UDC is G, implying Eq. 10 
with Nf = 2. The situation is more complicated, however, since 
we usually want the simulations to agree with the cosmological 
model of our observed universe, whose parameters are reasonably 
well known. Thus, some of these cosmological parameters might 
be treated as UDCs, depending on the purpose of the simulations. 

For example, if the Hubble constant Ho is treated as a UDC, 
then it would imply rj = 1, which together with G (that implies 
£ = a 3 ; Eq. (6)), results in A^ = 2, Nf = 1 and 

(m'o, l'o, t'o) = (a 3 m , al , t ) , 
Q' {t' = t,? = aR) = a B+3A Q(t, f) . (9) 

In particular this would leave the mass density unchanged, p'(f') = 
p'(t) = p(t), so that the effective f2 M (/) that is implied by the average 
value of (p) over the computational box would still follow the same 
original assumed cosmology. 

In cosmological A'-body simulations, the initial conditions are 
considered to be scale-invariant, since the amplitude of the initial 
fluctuations in the gravitational potential are (at least nearly) inde- 
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pendent of the wavenumber k. However, the corresponding fluctu- 
ations in density scale as (6p/p)(k) oc k 2 . This introduces a time 
dependent scale, li(t) = 2n/k\(f), at which the density fluctua- 
tions become of order unity, (8pl p)[k\(t)] = 1, and thus enter the 
strongly non-linear stage of their evolution. This scale changes un- 
der a rescaling of the length units, and so does <x 8 or the normaliza- 
tion of the initial power spectrum (i.e., the length scale that cr 8 rep- 
resents would no longer be 8/i~' Mpc upon rescaling of the length, 
but instead a x 8/7~' Mpc). Sticking to the observed value of cr 8 , 
and treating both Ho and <x 8 as UDCs, would effectively remove 
the last degree of freedom in our rescaling (i.e. result in Nf = 0). 
If, however, the cosmological parameters such as H or cr 8 are not 
treated as UDCs, and are allowed to vary (even if only over a lim- 
ited range that is consistent with current observational constraints), 
then a rescaling of the units given by Eq. l[7} could help to reduce 
the number of simulations required in order to numerically study a 
large parameter space with different cosmologies. 



3 CAVEATS 

Depending on the physics that are included in a simulation, further 
restrictions may arise in cases where there are additional UDCs or 
typical scales. In the previous section, the equation of state was 
implicitly assumed not to introduce any UDC, such as in the case 
of a simple polytropic equation of state, p = Kp y , where K can vary 
(with the specific entropy). However, this is not always the case. 

For example, degeneracy pressure in the Newtonian regime 
fixes the value of p/p 513 ~ % 2 /(m^m g) where m d (n d ) is the mass 
(number density) of the degenerate species while m eff = p/n d (here 
it is assumed that p is used as a primary hydrodynamic variable, 
rather than the number density n). This requires that ( = cfirf 3 . 
Since degeneracy pressure is usually important only when gravity 
plays a role as well, this would also require ( = cPrf 2 or altogether, 
r] = a 3 = f" 1 . In the relativistic regime, where the uncertainty 
principle implies relativistic velocities of the degenerate species, 
p/p* 13 ~ hcjin^ is fixed, implying £ = a?rf 6 . Together with grav- 
ity that introduces G, this implies ( = 1 and rj = a 312 . For ex- 
ample, the Chandrasekhar mass is approximately given by the 3/2 
power of the ratio of these two constants, Met, ~ (tic/m^G) 3/2 = 
Mp hnck /m^ ff . The transition between the two regimes of degeneracy 
pressure occurs when the mean distance between degenerate par- 
ticles is comparable to their Compton wavelength, thus fixing an 
absolute length-scale in the problem and requiring a = 1 if it ap- 
pears in the simulation. Together with gravity this would leave no 
free parameter (Nf = since N u d c = 3), and require ( = a = rj = 1. 

The optical depth t determines the probability for interaction, 
1 - e~ T , and must therefore remain unchanged. Thus, if p (rather 
than n) is a primary hydrodynamic variable then since dr = pK,dl, 
once the opacity coefficient k, of the matter is specified it should 
not change, and since it has units of P/m, this implies £ = a 2 . With 
the inclusion of radiation that introduces c as a UDC and requires 
a = 7], this implies f = a 2 = rj 2 . If, alternatively, n is the pri- 
mary hydrodynamic variable then since dr = cr,ndl this requires 
the cross-section cr, not to change and thus a = 1, which together 
with the inclusion of radiation (implying r] = a) gives or = r\ = 1 . 
If the mass of each particle is also to remain constant, this requires 
( = 1 leaving no degree of freedom and implying £ = a = r\ = 1. 

Optically thick radiation, or radiation pressure can also in- 
troduce UDCs. A black body, e.g., emits a power per unit area of 



Table 2. The dimensional-based scalings for the GRB afterglow syn- 
chrotron spectrum, in terms of ( and a = r\ (in column 3) or k and A (in col- 
umn 4), for the flux density within the different power-law segments (PLSs, 
Q — > F v ,a-^V,g; top part), the spectral break frequencies (Q — * vi-vn; 
middle part), and the flux density at the break frequencies (Q —> F Y> i—F V) u; 
bott om part). The notation for the different PLSs and break frequencies fol- 
low iGranoT&Saril J2002h . Column 2 gives the dependence of F v on v in 
each PLS for F V! \~F Vi q, and otherwise the relevant break frequencies. 
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Fy,3 
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^(2p+l)/2 Q ,(l-4p)/2 


K (p+2)/3 A (4p-l)/6 


Fy,4 
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£ 8+2/? a 8+2/? 
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Fyfi 
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5+2/? 5-2/? 
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Zp-5 
^2/3^30+6, 


Fy.l 


V7 


£2/5^4/5 


^-4/15 


fv,8 
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^■7/12^1/4 


^2/3^-1/12 


fv,9 
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f l/2 ffl l/2 


k 2 ' 3 A- 1 ' 6 


Fy,l0 


I'll) 


^-13/5^-18/5 


K V5tf/S 



ratio introduces c = 4<r/a (the radiation streaming velocity) that 
implies a = rj. Since k%T ss p ra d/« then p m d/(kgT) 4 that has unit 
of (ml 3 r 2 y 3 must also remain the same, implying f = rfa~ 3 , and 
together with the previous constraint, a = rj = f~ . If gravity is 
added as well then no freedom is left (Nf = and f = a = rj = 1). 

Radiation reaction (the force on accelerating charged parti- 
cles due to the back-reaction to the radiation they emit) or the ef- 
fects of radiative losses on the cooling of the radiating particles, 
can introduce additional dimensional parameters, that are universal 
for a given particle species, such as electrons, and thus introduce 
constraints on the scaling parameters. 



4 GRB AFTERGLOWS 

The dynamics of GRB jets during the afterglow stage have 
been numerically model ed using sp e cial r ela tivistic hydro- 
i Granot et aTT 1200 ll ; ICannizzo et al.l 



dynamic simulations 
|2004 IZhang & MacFadverJ 



2009, 2010; 



van Eerten et al 



err 4 and has a pressure of p ra( i = r\aT 4 , thus introducing the Stefan- 
Boltzmann constant cr = ac/4 and the radiation constant a. Their 



IWygoda. Waxman & Fraill 12011 



2009; Mimica, 



anniz zo et ; 
Giannios & Aloy 



Mimica, (jian mos & Aloy 
3; iMeliani & Keppensl 1201 ol ; 
De Colle et al. 2011a, 2011b). 



201 



As discussed above, for such simulations there is one UDC, c, 



© 201 1 RAS, MNRAS 000.[T]l6] 



Scaling relations in numerical simulations 5 



which implies a = rj and Eg. Q. It has recently been pointed ou0 
( Ivan Eerten. van der Horst & MacFadvenll201 ll) that the dynamics 
in this case obey a simple scaling relation, 



E' 
~E 



P -=A 



/' _ f _(K\ 1 ' 3 

1 ~ 7 ~ u; 



(10) 



which was justified by resorting to dimensionless or similarity vari- 
ables. However, this scaling simply arises from the freedom in the 
choice of the basic physical units, as described above. In particular, 
it corresponds to £ = k and a = rj = (k//1) 1/3 , or equivalently to 
k = £ and A = £/a 3 = '. This scaling holds regardless of the 
initial conditions or symmetry of the problem, and has nothing to 

do with self-similarity of the hydrodyna mics. 

It was also pointed out recently jyan Eerten & MacFadvenl 

l201lh that this scaling of the dynamic^ can also be extended to 
a similar scaling of the resulting afterglow synchrotron emission 
or the observed flux density F„, within each power-law segment 
(PLS) of the spectrum. This arises since within each PLS the local 
emissivity can be expressed as the product of a dimensional con- 
stant and a dimensionless function of the hydrodynamic quantities, 
so that a change in the basic units would affect only the dimen- 
sional constant, which would scale in a simple wayl 4 ] Therefore, 
such a rescaling of the basic physical units holds quite generally 
within each PLS, regardless of the dyna mics. In particular, the same 
rescaling holds in th e early relativ i stic dBlandfqrd & McKeelll976f) 
and late Newtonian JSedovlll946l :l Tavlorll 1950h (quasi-) spherical 
self-similar phases, as well as in the intermediate phase where the 
dynamics are not self-similar. Moreover, this scaling depends only 
on the PLS, and within a given PLS it does not depend on the exter- 
nal density profile (and would be the same for a uniform external 
medium and for a wind-like external medium). 

It has been demonstrated that when the dynamics are self- 
similar, a more elaborate scaling exists in which the flux density 
F v within each PLS s cales as a power-law with essentially all of 
the model parameters dSari. Piran & Narava n 1998; G ranot & S aril 
1 20021 ; Ivan Eerten & MacFadvenl 201 1 ). However, the dependences 
on the individual model param eters change between the relativis- 
tic UlandfordA^fei J1976I) and the Newtonian Sedov-Taylor 
self-similar regimes, and such simple power-law dependences on 
all of the model parameters do not exist in the intermediate phase, 
or whenever the dynamics are not self-similar. 

The dimensional-based scalings of F v hold only locally within 
each PLS, and change between different PLSs. This corresponds to 
different scalings for the break frequencies that separate the PLSs, 
so that their ratios changes under such a scaling, despite being di- 
mensionless. The lack of a global rescaling of the units for the ob- 
served radiation results in the need to parameterize and change "by 
hand" the spectral regime in order to calculate the lightcurve at 
a given observed frequency as it switches between different PLSs 
(when it is crossed by a break frequency). 

The lack of such a global rescaling of the units for F v can 
be understood as follows. Technically, it can be attributed to the 
fact that the local emissivity is separable, i.e. can be expressed as 
a product of a dimensional constant and a dimensionless function 
of the hydrodynamic variables, only within each PLS, and that this 



2 IScheck et allfepoA have outlined a similar scaling in a different context. 

3 There they use the scaling n' In = A for the number density, but this is 
effectively equivalent to p' /p = A since they assume that p/n = m p = const. 

4 According to the units of the part that scales with the hydrodynamics 
variables, and does not involve the distance from the source to the observer 
or UDCs such as the electron or proton mass or electric charge. 



dimensional constant that determines the scaling changes between 
different PLSs. The more basic reason behind this is that the emis- 
sion process introduces additional UDCs relative to the dynamics, 
even in the optically thin regime. For example, the radiation cares 
also about the total number of particles, i.e. about the number den- 
sity n and not only about the rest-mass density p, while p/n = m tS 
is usually taken to be constant (often set to the proton mass, m p ), 
and thus introduces a new UDC. Additional UDCs are introduced, 
e.g., through the synchrotron break frequencies, since they relate to 
the typical synchrotron frequency and cooling of the radiating rela- 
tivistic electrons, which have a universal mass and electric charge. 

Let us consider an emitting region of bulk Lorentz factor T, 
in the downstream region of a shock with a relative upstream to 
downstream Lorentz factor r u( j and upstream proper rest mass den- 
sity p u . For the afterglow forward shock p u is the external density 
and r u( i - 1 a T u j = T » 1 , while for the reverse shock p u is the 
density of the original outflow and typically T » r„d > r„d - 1 ~ 1. 
Thus, both shocks can be treated together. The comoving magnetic 
field scales as B 2 oc esr il( |(r u( j - l)p u , and the typical Lorentz fac- 
tor of the electron random motions scales as y„, ~ e e (m p /m e )(Tua - 
1) r u( j - 1. Thus, the typical synchrotron frequency scales as 



Y(eBlm e c)yf n oc e l J 2 e 2 TT 1 ^ (T uA - l) 5 ' 2 p u 



Hi, 



Pu 



1/2 



so that 



v' m /v m = ( >,2 a~ il2 — > A [/2 since y,„, T, T ua -, as well as the shock 
microphysics parameters e e and e B are all invariant under rescal- 
ings of the basic physical units that conserve c (a = rj). Note 
that the part involving UDCs, e/m e c, was not included in the scal- 
ing, since it is universal and does not change with the scaling of 
the hydrodynamic variables. The cooling break frequency scales as 



T(eBlm e c)y] oc I^JTV* oc r 1 [e s r ud (r ud - l)p u ] 



Pu 3/2? ob S ' and thuS V d v c = <T 



obs 

_ r-3/2 ~9/2„-2 



-3/2,-2 



r 3/2 a 5 ' 2 -> k- 2 ' 3 A- 5 ' 6 , 
where r obs is the observed time (when the emitted photons reach the 
observer) and y c = 6mn e c/(cr T B 2 Tt obs ) oc B~ 2 r _1 r? is the random 
Lorentz factor to which the electrons cool on the dynamical time. 
Again, parts involving UDCs, such as e/m e c or m e c/cr T , were not 
included in the scaling. A global rescaling of the units would re- 
quire v m p„ 1/2 and v c pu /2 r 2 bs with units of m~ ll2 P l2 r 1 and m 3/2 r 9/2 t, 



i = a 3 rj- 2 ' 3 



v 7 '- 5 



, or altogether <f = a = rj = 1, which eliminates 
all of the freedom in such a rescaling. The peak synchrotron flux 
density scales as F„_ raax oc TBN e oc e^rr^ 2 ^ - l) 1/2 p u /2 M oc 
p'J 2 M (where N e and M are, respectively, the isotropic equivalent 
number of emitting electrons and rest mass in the shocked region, 
and M/N e = m eff = const), which implies that -FJ, max /F vmilx = 
£3/2 ff -3/2 _^ K ^}l 2 , Note that the distance to the observer, D, is not 
included in the scaling of F vmax since it does not change with the 
hydrodynamic variables. 

Even though there is no non-trivial global scaling of the units 
that obeys Eq. {3j f° r the flux density (Q — > F v ), such a scaling still 
works locally within each PLS (labeled by a subscript T), 



KjOL = ttf obs) = i ai a b 'F VJ (t obs ) , 



(11) 



where the dependence on t and r is replaced by / obs . This can be 
understood since the flux density within each PLS is the prod- 
uct of F vmax and certain fixed powers of the break frequencies 
(v m , v c , and the self-absorption frequency that is not discussed 
here for simplicity), whose scalings can be derived from sim- 
ple dimensional considerations (as shown above). For example, 
F v ,d ~ F„ raax (v/v„,) l/3 and F„, F x F v , max (v/v c )~ 1/2 for PLSs D 
and F, respectively, using the notations of iGranot & Saril d2002l) . 
This implies that F' vD l F vD = f^ar 1 ^ kA 1/3 (a D = 4/3 and 
b D = -1) and F' vF /F,, F = f 3/4 o<- 1/4 -> K 2j3 A lln (a F = 3/4 
and b F = -1/4). As illustrative examples of how the scalings for 
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self-absorbed PLSs may be derived, one can readily obtain that 
F„, B x niRjYDfilv 2 lc 2 )Ty m m e c 2 oc v 2 R 2 implying F' yB /F v , B = 
{°a 2 -> k 2/3 A- 2 > 3 (a B = and b B = 2), while for PLS A 
y m is replaced by y e (y) oc (v/TB) >/2 [obtained from requiring 
v ~ v syn (y e ) ~ T(eB/m e c)y 2 ], implying F vA oc y 5 ' 2 ,/^" 4 and 
KJ F -«A = r 1/4 a 11/4 -» K 2/3 A- n,i2 (a A = -1/4 and b A = 11/4). 
Therefore, these scalings (or a, and b,) do not depend on the ex- 
ternal density profile or on the details of the dynamics (and are the 
same in the relativistic and Newtonian self-similar regimes, when 
the dynamics are not self-similar, or for the reverse shock). All of 
the different scalings are summarized in table [2] 



5 DISCUSSION 

The freedom in the choice of units in the dynamical equations that 
describe the evolution of different types of physical systems and in 
their solutions, has been outlined and elucidated. The main results 
are summarized in Table [TJ While the emphasis was on numerical 
solutions of the dynamical equations through simulations, similar 
scalings hold equally well for analytic solutions of the same equa- 
tions. The number of free parameters Nf that describe the family 
of physical systems that corresponds to a given solution of such a 
set of equations is given by max(0, 3 - iV udc ) (Eq. |2)), where iV u( jc 
is the number of independent (in terms of their units) universal di- 
mensional constants (UDCs, such as c, G, h, m e , etc.). This corre- 
sponds to the three basic physical units (of mass, length and time) 
while accounting for the independent constraints on their possible 
rescalings. Such rescalings of the basic units are potentially rel- 
evant to many different areas of research, such as plasma physics, 
astrophysics, cosmology, fluid dynamics or Earth and planetary sci- 
ences. They can prove very useful in numerical studies of various 
physical systems, and save precious computational resources, espe- 
cially in systematic numerical studies of a large parameter space. 

The author thanks E. Nakar, F. van den Bosch, F. De Colle, T. 
Piran, O. Bromberg, E. Ramirez-Ruiz and the anonymous referee 
for useful discussions, suggestions or comments on the manuscript. 
This research was supported by the ERC advanced research grant 
"GRBs". 
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